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Abstract. We present an introduction to the backflow effect in quantum mechanics - the 
phenomenon in which a state consisting entirely of positive momenta may have negative current 
and the probability flows in the opposite direction to the momentum. We show that the 
effect is present even for simple states consisting of superpositions of gaussian wave packets, 
although the size of the effect is small. Inspired by the numerical results of Penz et al, we 
present a wave function whose current at any time may be computed analytically and which 
has periods of significant backflow, with a backwards flux equal to about 70 percent of the 
maximum possible backflow, a dimensionless number ctm ~ 0.04, discovered by Bracken and 
Melloy. This number has the unusual property of being independent of h (and also of all 
other parameters of the model), despite corresponding to a quantum-mechanical effect, and we 
shed some light on this surprising property by considering the classical limit of backflow. We 
conclude by discussing a specific measurement model in which backflow may be identified in 
certain measurable probabilities. 



1. Introduction 

A striking but little- appreciated phenomenon in quantum mechanics is the backflow effect. This 
is the fact that, for a free particle described by a wave function centred in x < consisting 
entirely of positive momenta, the probability of remaining in x < may nevertheless increase 
with time. That is, the quantum-mechanical current at the origin can be negative and the 
probability flows "backwards". 

This effect was first noted by Allcock in his early work on the arrival time problem in quantum 
theory [IJ, but was not studied in detail until 1994 when Bracken and Melloy carried out the first 
systematic study [2j. Their most important discovery was that there exists a limit on the total 
amount of backfiow. Backfiow means that the probability of remaining in x < may increase 
with time, but this increase is bounded by an amount Cbm, a dimensionless number computed 
numerically to be approximately 0.04. Furthermore although backflow is clearly a nonclassical 
effect, the quantity Cbm is independent of h (and also of the mass of the particle and the time 
duration of the effect.) For this reason Cbm has been described as a new quantum number. 

The effect has been further investigated by a number of authors. Better estimates for Cbm 
were given in Refs.^ [4], with the latter also giving a numerical estimate of the maximum 
backflow state. Some analytic examples of backflow states and the measurability of the effect 
were explored in Ref . |5j . Connections between backfiow and the phenomenon of superoscillations 
were noted in Ref. [6j. The effect has also been explored in the Dirac equation for a particle 
in a linear potential [8] and for angular momentum [9]. Most recently the present authors have 



investigated a number of aspects of the backflow effect, including the classical limit [10]. Finally, 
the effect is often noted in connection with the arrival time problem in quantum theory 

Apart from the studies cited above, there has been little discussion of the effect in the 
literature to date. The purpose of this contribution is to provide an introduction to the backfow 
effect both as a reference for researchers encountering it in their studies of quantum theory and 
also hopefully as an invitation to encourage others to take an interest in the effect. 

We begin in Section 2 with a detailed formulation of the problem. We define the current and 
the flux and consider the properties of the spectrum of the flux operator, which allows us to give 
the clearest statement of the backflow effect. 

In Section 3 we give a simple example of a state with backflow consisting of a superposition 
of gaussian states. The backflow for such a state is, however, rather small. 

In Section 4, we review the numerical computation of the maximal backflow state and 
eigenvalue. We also present an analytic wave function with backflow which appears to match 
closely the numerical solutions for the maximal backflow state by Penz et al [4J. The current 
at arbitrary times for this wave function may be computed analytically and we find that it has 
a backflow of approximately 70 percent of the maximal value. This is a much larger backflow 
than any analytically tractable states previously discovered. 

In Section 5, we consider the naive classical limit 7i — )• of backflow, and in particular, 
we address the fact that the bound on backflow Chm discovered by Bracken and Melloy 
appears to be independent of h. We show that the expected dependency on h reappears in 
realistic measurement models, where measurements are described not by exact projectors but 
by quasiprojectors involving parameters characterizing the imprecision of real measurements. 
Under these conditions the naive classical limit is restored. 

In Section 6, we consider a simple measurement model and discuss the ways in which backflow 
may be seen in the probabilities for measurements. 

We summarize and conclude in Section 7. 

2. Detailed Formulation of the Problem 

2.1. The Flux 

We consider a free particle with initial wave function 'ip{x) concentrated in x < and consisting 
entirely of positive momenta. We consider the probability flux F{ti,t2) crossing the origin during 
the time interval [ti,t2], defined by 




(1) 



where J{t) is the usual quantum-mechanical current at the origin 
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dx dx 
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The fiux may also be written in terms of the Wigner function [12j at time t, Wt{p, q) 




(3) 



(For a useful review of the properties of the current and related phase space distribution 
functions, see Ref.[T3j). It is useful to write these expressions in an operator form. We introduce 



projection operators onto the positive and negative x-axis, P = 6{x), and, P = 1 — P = 9{—x) 
respectively. The flux may then be written in terms of the operator F{ti,t2) deflned by 



F{h,t2) 



P{t2)-P{tl) 
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where the current operator is given by 



J = — {'p5{x) + 5{x)p) 



(5) 



So Eq.([T]) may also be written 



F{ti,t2) = {Fitut2)) = [ 



dt (V|J(t)|^>. 



(6) 



The flux Eq.([T]) is a difference between two probabilities and is positive when those probabilities 
behave according to classical intuition, i.e., when the probability of remaining in x < decreases 
monotonically. 

As indicated, in the full quantum-mechanical case, the flux can be negative. One way to see 
this is to note that the Wigner function need not be positive for general states |12) . and thus 
Eq.(l3|) need not be either. Another way is from the current operator Eq.([5]): both p and 6{x) 
are non-negative operators on states with positive momentum, but since they do not commute, 
the current operator J is not a positive operator. 

2.2. Most Negative Flux as an Eigenvalue Problem 

Following Bracken and Melloy [2], a useful way to investigate the backflow effect is to look at 
the spectrum of the flux operator Eq.dH) (restricted to positive momenta). We thus look for 
solutions to the eigenvalue problem 



where the states |$) consist only of positive momenta. (We choose an opposite sign convention 
to [2] so that the backflow states have A < 0). The most negative value of the flux F{ti,t2) is 
then given by the most negative eigenvalue of Eq.Q. 

By time evolving the state we may choose the time interval [ti, to be [— T/2, T/2] and the 
eigenvalue equation in momentum space then reads 



eip)Fitut2m = xm 



(7) 




Defining rescaled variables u and v hy p = 2y^mh/Tu and k 



1^ mh/Tv Eq.® becomes 




(9) 



where 4'{u) = {mTi/ATY/^<^{p) and is dimensionless. Crucially all physical constants have 
dropped out of this equation so that the A are dimensionless and independent of fi, m and T. 



This eigenvalue equation has been studied by a number of authors, both numerically and 
analytically [21 ElS] . The eigenvalues lie in the range 



- Cb„ < A < 1 (10) 

where was computed numerically and found to be 

Cbm 0.038452. (11) 

It was conjectured in Ref.[2j that the spectrum is discrete in the interval [—Cbm, 0] but continuous 
in the interval [0, 1]. The extremizing state was given numerically by Penz et al [1]. We will 
display this maximising state and also a good analytic expression approximating the numerical 
results in what follows. At present there is no analytic account of the properties Eqs. ()10p . ()lip . 
Note that the eigenvalues are independent of T which means that the duration of a period of 
backflow can be arbitrarily long. 

3. Backflovif for superpositions of gaussians 

In this section we show that the backflow effect arises in the familiar, and also potentially 
experimentally realisable, setting of a superposition of two gaussian wavepackets. However 
gaussian wavepackets have support on both positive and negative momentum, so we will also 
have to show that this is not the result of any initially negative momentum. We will see that 
whilst superpositions of gaussian states do indeed give rise to backflow, the size of the effect is 
considerably smaller than the theoretical maximum. In this Section we work in units in which 
h = 1 and we set the particle mass ra = 1. 

Bracken and Melloy demonstrated that the backflow effect may be observed in the simple 
case of a superposition of two plane waves [2]. This state is unnormalisatble, but it can be 
turned into a more physical state by replacing the plane waves with gaussians tightly peaked 
in momentum, without affecting the basic conclusion that the state displays backflow for well 
chosen values of the various parameters. 

Consider the normalised state 

V^(x,i) = E^A,-^=i=exp {^p,{x-p,t) - ^^^^) ■ (12) 

This is a sum of two initial gaussian wavepackets with equal spatial width cr, evolved for a time 
t. If we let (T — 7- oo we essentially recover a sum of two plane waves. Instead of giving the 
rather complex expression for the current, we will simply plot the current at the origin and 
the probability of remaining in x < as functions of time for the state in Eq. (|12|) and for the 
following set of parameters; 

pi = 0.3, P2 = 1.4, a = 10, Al = 1.8, A2 = 1. (13) 

We clearly see from these plots that there are several intervals during which the current is 
negative. A magnification of one of these backflow region is shown in Fig. ([3]). The effect is 
robust with respect to small changes to these parameters. 

The parameters in Eq. ()13p give rise to the greatest amount of backflow we have been able to 
find, although we have not searched the entire parameter space. The value of the flux during 
the largest period of backflow is 





Figure 1. Plot of the current for a 
wavefunction consisting of a superposition of 
two gaussians, with the parameters given in 
Eq.dn]). 



Figure 2. Plot of the probability for 
remaining in < for a wavefunction 
consisting of a superposition of two gaussians, 
with the parameters given in Eq. (jl3p . 
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Figure 3. Close up of Fig.©. P{t) is clearly 
seen to increase between t « 2 and t ^ A. 



where J{t) is the current, Eq.(l2|), and the interval [ti,t2] is chosen such that the current is 
negative for the whole of this time. Computing the flux during this interval gives, 

-0.0061, (15) 

or about 16 percent of Cbm- 

It is important to check that this probability backflow cannot be explained by the tiny 
probability of having negative momentum in this gaussian state. It can be shown that the 
probability that a measurement of the momentum of this state would yield a negative answer is 
of order 10"^", so the negative flow of probability is entirely due to the backflow effect [lOj. 

4. An approximation to the backflow maximising state 

Backflow states may be found by solving the eigenvalue equation Eq.Q. The numerical work 
of Penz et al [1] yielded an approximate eigenstate satisfying Eq.Q giving the most negative 
eigenvalue —Cbm, i-e., the largest amount of backflow. It is of interest to find analytic expressions 
for wave functions matching these numerical results closely and for which the current may be 
computed analytically. We give such a state in this section. To be clear, we are not giving 



an approximate analytic solution to the eigenvalue problem, Eq.(l9]). Rather, inspired by the 
numerical solution, we will exhibit an analytic expression for a wave function which closely 
matches the backflow maximising state, and show that is has significant negative flux. 



4-1. Numerical results 

We first review the numerical results of the computation of the backflow state. We have repeated 
the numerical analysis of Penz et al.[¥] of the optimizing state, (j)max, and its current in order 
to compare with our analytic approximation. In Fig.dU we plot the numerically computed 
maximum backflow state, (t)max, together with the function 

, , , 1 cos(u^) , , 

<t^as{u) = (16) 

which seems to match well the asymptotic form of cjjmax ■ 




We plot in Fig. ([5]) the current computed from (pmax- Note that the current appears to have 
a specific singularity structure at t = ±1, where it jumps between ±00. This is presumably 
related to some properties of the flux operator, but it also seems to depend on the asymptotic 
behavior of (jimax [13] • These two plots are in agreement, in general shape, with the numerical 
results of Penz et al [3] and we will compare our analytic results with these plots in what follows. 



4-2. Analytic Approximation 

Consider the momentum space wavefunction 



ae-^ + Q- C{u)) 
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(17) 



where is a normalisation factor. Here 



C(u) = FresnelC 
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2 r 



dx cos(2;^). 



(18) 



This has the asymptotic form. 



4>Aiu) ~ N 
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(19) 



Note that this does not in fact match the asymptotic form of 
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Figure 6. (I)a{u), for a = 0.6, b = 2.8 (solid Figure 7. Current, J{t) for a = 0.6 and 
line), with cj)max{u) for comparison (dashed 6 = 2.8. 
line) . 

We plot (j)A in Fig. ([6]) for the values of a and b which maximize backflow. We see good 
agreement with the numerical result. The maximum amount of negative flux we can generate is 



which occurs for the parameters a = 0.6, b = 2.8. This corresponds to about 70 percent of 
Cbm- We plot the current J{t) for these parameters in Fig.Q. The current is close to the 
numerical result. Fig. ([5]) away from t = ±1 but lacks the correct behavior as t — ?■ ±1. This 
is presumably related to the fact that the asymptotic behaviour of (pA does not match that of 
(pmax- It would be interesting to know if there exists a better analytic approximation to the 
maximum backflow state that does have the correct asymptotic behavior and still gives rise to 
an analytically computable current. This will be explored elsewhere [14] . 

5. The Classical Limit of Backflow 

Some insights into the properties of backflow may be found by looking at its classical limit. 
Strictly one should do this via the usual framework of open system dynamics, and a detailed 
discussion of the current for an open system was given in Ref.|15j. In this work the resulting 
positivity of the current, after finite time, is clearly seen. 

However it is nevertheless interesting to look at the naive classical limit of backflow 
characterized by the limit ^ — )• 0. As seen from Eq.(l9]) the eigenvalues of the flux operator 
are independent of h. Thus although the existence of negative eigenvalues (negative flux) is 
clearly a quantum phenomenon, backflow does not go away in the naive classical limit, ^ — )• 0, 
as one might have expected. The origin of the independence of the eigenvalues of the flux 
operator on fi is the fact that there is no way to construct a dimensionless number from the 
parameters in the problem, h, m and T. What is missing is a length scale. 

In any realistic experimental set up, measurements have finite resolution. It is therefore often 
more appropriate to model the measurement process using POVMs rather than exact projection 
operators. With this in mind, instead of defining the flux operator in terms of exact projection 
operators P = 9{x), instead let us define it in terms of a quasiprojector Q. This seems reasonable 
since, as discussed earlier, backflow can be measured by measuring whether the particle is in 
a; > at two different times and, due to the inevitable imprecision of real measurements, such 
measurements are best modeled by quasiprojectors. A convenient choice of quasiprojector is 



F = - 



0.02757, 



(20) 




(21) 



where 5a{x — y) is a smoothed out (5-function, 
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(22) 



This goes to the usual (5-function as u — t- and then Q — )• 9{x). By replacing the exact projector 
in the definition of the fiux with this new quasi-projector we may compute modified expressions 
for the current and the flux |10] . In particular the eigenvalue equation Eq.([9]) becomes, 



1 



vr 



dv 



sin(u2 — t>2) 
(u — v) 



(23) 



where the dimensionless number a is given by a2 = 2ma'^/hT. This means that the eigenvalues 
A will now depend on a, so we write A = A(a). Through a they will therefore depend on h and 
the "limit" 7i — ?• now clearly means the regime a ^ 1, that is, h <C 2mcj2/T. Hence, in a more 
realistic measurement situation, the bound on the total backflow - the most negative eigenvalue 
of Eq. ()23p - will depend on h and the limit ^ — )• may now be more meaningful. 

A reasonable conjecture is that the negative eigenvalues will increase with a and also that 



A(a) > -Cbr, 



(24) 



for all a, so that —Cbm emerges as a lower bound on the eigenvalues, achievable only in the limit 
a — )• 0. It seems unlikely, however, that all the negative eigenvalues will all become positive or 
zero, except perhaps in the limit a — )• oo. We have not been able to solve Eq. ()23p analytically, 
so instead we have obtained numerical estimates for A(a) for various values of a, and we plot the 
result in Fig.®. The value of A (a) does indeed increase with o, tending to zero asymptotically. 
In fact, numerical solutions are consistent with the asymptotic form. 
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A(a) 



(25) 



for large a, which can be understood by examining the asymptotic form of Eq. (j23p for large a 

m- 
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Figure 8. The most negative eigenvalue of 
Eq. ([25|) . A(a) as a function of a. 



In summary, for quasi-projectors, Eq. ()2ip . which are more realistic models for measurements 
than exact projectors, the eigenvalues, in particular the lowest eigenvalue representing the most 
negative flux, do depend on h (and m and T). The lowest eigenvalue appears to go to zero for 
a — )• oo. This indicates that all negative eigenvalues go to zero (or become positive) in the naive 
classical limit /i — )• 0. By contrast the positive eigenvalues are not significantly affected [lOj . 



6. Backflow and Measurement Models 

In this section we relate the above results on backflow to measurements. The aim is to begin to 
address the practical question of how backflow may be measured. 

Many more elaborate and realistic models for the measurement of the arrival time (involving 
model detectors, for example) naturally lead to an arrival time probability deflned with a complex 
potential. This is described in detail in many places |16) . These models typically yield an arrival 
time probability distribution which is closely related to the current and from which the current 
may be extracted, even when negative, thereby leading to a possible measurement of backflow. 

Consider an initial wave packet starting in x < with positive momentum. We seek the 
arrival time probability distribution Il{T)dT for crossing the origin between r and r + dr. We 
consider a complex absorbing potential of step function form in x > so the Hamiltonian is 
Hq — iVQ9{x), where Hq is the free Hamiltonian. We define the survival probability N(t) to be 
the norm of the state at time r after evolution with this complex Hamiltonian. Then 



This expression may be simplified considerably in the usual weak measurement approximation 
(Vo small compared with the energy of the state), yielding |10j . 



where \ipt) = e~*^°*|^). This is the expected semiclassical result [I1|T71[T6]. Note that Eq. ([27|l 
is not necessarily positive, due to the negativity of the current in certain states. However this 
is an artifact of the approximation and should not matter for sufficiently small Vq, and we will 
assume that Eq. ([27|) is positive. 

The quantity n(r) corresponds to the arrival time distribution measured by a realistic 
measurement so can in principle be determined experimentally. The current can then be 
extracted from Eq. ()27p by deconvolution [17]. This therefore gives a method of measuring the 
current and the flux, and checking for backflow. In Fig. ([9]) we plot the measurement probability 
Eq. ()27p for two values of Vq and also the original numerically computed current, to see how 
the time-smearing affects the backflow. We see that positive regions of the current are not 
qualitatively changed very much, in keeping with semiclassical expectations, but negative regions 
of the current become positive as a result of the smearing, as they must, since the measured 
probability is positive. Details of other possible measurements of backflow are given in Ref.|lUj. 

7. Summary and Conclusions 

The purpose of this paper was to introduce various aspects of the backflow effect. After setting 
up the problem in Section 2, in Section 3, we showed that backflow may occur for a superposition 
of gaussian wavepackets. These states are important since they are experimentally realizable, 
but the maximum amount of flux is very small, only about 16 percent of the maximum possible. 

In Section 4, we gave an analytic expression for a state matching rather closely the numerically 
computed states giving maximal backflow, computed by Penz et al [5]. The plot of the current 
shows reasonably good agreement with the numerical solution, except at the end points t = ±1 
of the backflow region. We computed the most negative flux for this state and found it to be 
about 70 percent of the numerically computed maximum backflow, significantly better than any 
previous analytic expression for a backfiow state. 

In Section 5, we discussed the classical limit of backfiow. The most interesting aspect of this 
is the issue is that the eigenvalues of the flux operator are independent of h. This appears to 



n(r) 



dN 




(26) 




(27) 




Figure 9. A plot of the current (solid line) 
and time-smeared current Eq. ()27p for Vq = 
0.5 (dashed line) and Vq = 0.1 (dotted line). 



mean that backflow does not go away in the naive classical limit /i — t- 0. We showed that this 
situation starts to appear more reasonable when the projectors used in the definition of the flux 
operator are replaced by quasiprojectors, which include a physical parameter characterizing the 
imprecision of real measurements. The eigenvalues then do depend on h and there is evidence 
that all the negative eigenvalues become zero or positive as ?i — t- 0, restoring the naive classical 
limit. However, there are clearly more issues to explore around this question. 

In Section 6 we discussed a simple measurement model in which the current can be obtained 
from the measured probability by deconvolution, and from this result the negative current could 
in principle be obtained. The features of backflow elucidated here may be of value in designing 
experiments to test backflow. These and related ideas with be explored elsewhere. 
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